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1. Introduction 

String theory outside the critical dimension has proven a fruitful subject of investiga- 
tion, yielding many valuable insights. The introduction of such "non-critical" strings by 
Polyakov [|l] highlighted the importance of local scale invariance as a consistency condition 
on the worldsheet QFT describing string propagation. Indeed, noncritical string models 
are now understood as a particular class of backgrounds of ordinary, "critical" string theory 
with asymptotically linear dilaton: 

G^^f]^ , $~n-X (1.1) 

where n 2 = 2 ^.~P for the bosonic string, and n 2 = 10 ,~P for the fermionic string. Per- 
turbative string theory in such a background is self-consistent if the dynamics avoids the 
region of strong coupling. One way to accomplish this (for D < 2) is to turn on a closed 
string tachyon background 

T = {ie an - X . (1.2) 

The resulting worldsheet potential pushes strings away from strong coupling. 

The formulation of non-critical string models in target dimension D < 2 in terms of in- 
tegrals over random matrices [@,|3],§] provided the first insights into stringy non-perturbative 
effects HH . Indeed, recent work [|7|,|8|,p|, |T0| , p!T| , |T2| has established that the matrix model de- 
scribes an open string dual of non-critical string theory background (|l.l| )-( [L2j ), in which 
the matrix eigenvalues are D-branes of the non-critical string. 

Thus these string backgrounds provide particularly tractable examples of the corre- 
spondence between gauge theory and gravity uncovered in recent years JT3J (for a review, 
see WW)- Here, both sides of the correspondence are amenable to precise calculation. 
Namely, on the matrix model/open string/gauge theory side, the matrix path integral 
is exactly solvable. On the closed string/gravity side, the correlation functions of the 
worldsheet CFT of the non-critical string can be computed exactly using conformal boot- 
strap methods []T5|JT^ , p!7| , p!8| , p!9| (for a review, see ||20|| ). The agreement between the two 



approaches of all quantities thus far computed provides strong evidence that the two for- 
mulations are equivalent. 

On the closed string side, another interesting non-critical string background is the 2D 



black hole 



k 

ds 2 = — \dr 2 - tanh 2 r dt 2 ] 

a' 1 J (1.3) 

$ = $ - 2 log [cosh r] 
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with k = 9/4 for the bosonic string, and k = 5/2 for the fermionic string. This geometry is 
in the class of asymptotically linear dilaton backgrounds , and the Euclidean geometry 
admits an exact worldsheet CFT description as the gauged WZW model SL(2, M)/C/(l). 
These results raise the possibility of a non-perturbative and perhaps solvable matrix model 
formulation of 2D black hole dynamics. 

An obstacle to the realization of this idea in the bosonic string is the non-perturbative 
ambiguity in the definition of the corresponding matrix model. In the latter, the matrix 
eigenvalues are free fermions in an inverted quadratic potential, with only the metastable 
states on one side of the barrier filled (see figure la). The interpretation of eigenvalues 
in the "ground state" tunneling to the other side of the barrier, or simply vaulting over 
the barrier in a high energy process as in figure lb, remained unclear. Since black hole 
formation is expected to arise through such high energy processes (c.f. |23j for a discussion 
in the present context), the status of black holes in the matrix model also remained unclear. 





Figure 1. The matrix model for the bosonic string: (a) The metastable "ground 
state"; (b) The puzzle introduced by high energy scattering. 

Now it is understood that the stable ground state, where eigenvalue fermions fill both 
sides of the harmonic barrier, is a matrix model formulation of the superstring - specifically, 
its type OB incarnation |TT|,OJ (see figure 2a). It seems appropriate to revisit the question 



of black hole formation in the matrix model. 




type OA 



Figure 2. The matrix model eigenvalue potentials and ground state configurations 
for the (a) type OB, and (b) type OA fermionic string. 



For reasons of technical simplification, we will consider the related type OA matrix 
model, which allows us to sidestep various subtleties associated to the intercommunication 
of the two sides of the harmonic barrier in the OB model The eigenvalue potential of 
the OA model 



V 



X 2 



+ 



'M) 



4a' ' 2A 2 ' 

depicted in figure 2b, has eigenvalues restricted to the half-line A < 0. The parameter q is 
the RR charge of the background. For q = the black hole solution of (|1.3|) applies (with 
k = 5/2), and can again be described on the worldsheet as the Af = 1 supersymmetric 
SL(2, JR)/U(1) coset model. The more general black holes with g ^ have been studied 
in m,E^]2|]. 

Can one form black holes in 2D string theory? In any situation where black holes form, 
one expects them to constitute the generic intermediate states due to their large density of 
states, and therefore to dominate the behavior of the S-matrix. Thus we would like to ask 
whether the S-matrix of 2D string theory has any characteristics that one would associate 
with the appearance of black holes as intermediate states. In higher dimensions, this is 
indeed the case in semiclassical gravity; the generic result of a very energetic collision is 
the formation of a black hole. 

Nevertheless, two arguments suggest that black holes do not form in 2D string theory. 
1. First, the matrix model has infinitely many conserved quantities, since the underlying 
dynamical objects are free fermions. Any black hole will correspond to rather specific 
values of these conserved quantities, and therefore the generic initial state having 
other values for these quantities will have small overlap with the black hole. 
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2. In other examples of the gauge/gravity correspondence, the formation of black holes is 
associated with a deconfinement transition in the gauge theory, in which one accesses 
non-singlet degrees of freedom in the gauge theory [^0[ . These degrees of freedom are 
projected out of the matrix model for 2D string theory, suggesting that there is no 
regime in which the dominant configurations in the dynamics can be associated to 
black holes. 

These two issues are related. The appearance of free fermions and their associated hierarchy 
of conserved quantities is directly related to the absence of non-singlet degrees of freedom 
in the matrix model, which would be needed to provide the exponentially large degeneracy 
of black hole states. The projection onto singlets is consistent with the Gauss law of the 
gauge theory on DO-branes, which provide the worldsheet description of matrix eigenvalues 
pulled out of the Fermi sea [12]. 



Even so, it would be nice to verify this conventional wisdom with a concrete calcu- 
lation, to see if there are any indications of an object with the characteristics of a black 
hole which participates as an intermediate state in the dynamics. For example, it might be 
possible for a forming black hole to shed via tachyon radiation the collection of conserved 
charges along with some fraction of the infall energy, leaving behind a long-lived black hole 
intermediate state. 

Below we will examine high energy scattering processes in the matrix model, using 
the formalism of for the exact S-matrix. Of course, the high energy behavior of the 
S-matrix is of intrinsic interest, independent of the question of whether black holes are 
involved in the scattering. We wish to look for the probability to find N particles in the 
final state, to see if there is a correspondence with the expected decay spectrum of a 2D 
black hole. We will examine in detail the 1 — > N scattering amplitude, but we expect other 
amplitudes involving a finite number of high energy incoming particles to behave similarly. 
Is there any feature in the outgoing spectrum that we might identify with a black hole 
intermediate state? We will show that, for N more than a few, the probability to produce 
iV particles is a rapidly decreasing function of N, with no feature in the range iV ~ ui that 
would be characteristic of black hole intermediate states. 



1 Even if non-singlet states were somehow allowed, their energy cost becomes infinite in the 



continuum limit of the matrix model [31|. 



2. The black hole in type OA gravity 



Before embarking on an analysis of scattering in the matrix model, let us consider what 
result would be predicted by the appearance of 2D black holes as the generic intermediate 
state. 

At lowest order in a', the effective action of type OA string theory admits a class of RR 
charged black hole solutions [26,^7,^8,2^]. In Schwarzchild-like coordinates, the geometry 
can be written (for \x = 0, i.e. no tachyon field background) 



ds 2 = -l(cj))dt 2 + 



$ = KC 



(2.1) 



where k = J — , and l(<fi) is given by 



1(4)) = 1 + e 2 ^-^) \2 K e- s ((j) - H ) - 1 



(2.2) 



Here, (p = </> H is the location of the black hole horizon; the horizon of the extremal black 
hole is located at 



1 

2k 



log 



32tt 
q 2 



and 



e = 2k( 



(2.3) 
(2.4) 



is a measure of the departure from extremality. 

In terms of these quantities, the Hawking temperature and ADM mass above ex- 
tremality are given by 



K 



T=—(l-e- s ) 

Z7T 



-E'adm — 2KC 



-2K6 



e s - 1 - e 



(2.5) 



Note that the temperature rises smoothly from zero at extremality to a limiting tempera- 
ture of order one in string units. The high mass behavior is thus Hagedorn thermodynam- 
ics, with an entropy 

S BH ~ (3 H M (2.6) 



where (3 H = 2-K^a' j2. 

The Hawking temperature is related to the actual spectrum of radiation seen at infinity 
via so-called greybody factors (c.f. [|33| ). The potential seen by a Hawking quantum as it 



travels away from the horizon acts as a filter and distorts the spectrum. Let us estimate 
the magnitude of such effects. 

The equation of motion for the closed string tachyon T = e~®T in this background is 



l^- ) T+{u' z -V)T = Q 



(2.7) 



where 



V = lull' - 15n 2 l(l - 1) 



2, 2k(0-0 h ) 



k le 



1 + 14e" £ - 2ne~ £ (0 - <£ H ) 



(2- 



A sketch of this effective potential is shown in figure 3. 



Figure 3. The effective potential seen by a tachyon perturbation in the OA black 
hole background, as a function of the "tortoise coordinate" x = f l~ 1 d(j). 

The height of the effective potential is of order one in string units. The transmission 
amplitude rises to unity for incident energies far above the barrier height, and vanishes 
smoothly as the incident energy goes to zero. At any rate, for black holes of mass above 
extremality E ADM ^> k, where the Hawking temperature is of order T ~ k/2-k, the greybody 
effects amount to a factor of order one. It is only for near-extremal black holes, where 
the typical Hawking quantum has an energy small compared to the barrier height, that 
greybody effects will be significant. 

One can also consider black hole formation in the context of the low-energy effective 
field theory. The gravitational back-reaction of an infalling tachyon pulse propagating 
according to (|2.7| ) naively results in black hole formation for any sufficiently energetic 
pulse.i The horizon moves out to large cf>, and therefore weak coupling, with increasing 



2 As in the study of dilaton gravity coupled to conformal matter [ 34 1 , we would have to consider 
a incoming tachyon pulse localized along Z~ in order to have a finite energy density which is 
sufficiently localized to create a black hole. This can be achieved by making a tachyon wave 
packet instead of the plane wave state which we will consider in the following sections. However, 
this does not change the qualitative feature of the number distribution in the final state discussed 
in section 5 (see also footnote 4). We would like to thank A. Strominger for raising this issue. 
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energy. The bulk of the pulse is not backscattered by the effective potential V, equation 
( |2.8|) . Provided that a' corrections to the low-energy field equations do not change the 
qualitative structure deduced from the leading order equations, the black hole is the generic 
intermediate state for any highly energetic initial state - once tachyons fall inside a trapped 
surface, the black hole will form. The trapped surface forms at weak coupling for sufficiently 
large energy, and so the classical field equations are reliable. 

The qualitative picture that then emerges is as follows: Energetic incoming states 
form black holes. For black hole masses far above extremality, E ADM ^> nq 2 , the Hawking 
temperature T H is of order one in string units. The decay of the black hole will proceed 
smoothly via the emission of quanta whose energy is of order one in string units, until 
E A dm ~ KQ 2 - At this point the character of the decay changes; the temperature varies 
with the mass, smoothly decreasing with as the black hole evaporates toward extremality. 
We can make these endpoint effects parametrically small by considering black holes whose 
initial mass is large compared with the charge. The expectation then is that the typical 
number of quanta emitted in a black hole decay will be of order the mass in string units, 

W^ = %^. (2.9) 

-Z H K/ZTT 

It should however be emphasized that the geometry ( |2.1| )-( pT2| ) is strongly curved in 
string units, and the above picture of its properties could be strongly corrected by stringy 
effects. It is somewhat reassuring in this regard, that for q = the black hole admits an 
exact CFT description as an SL(2, H)/C/(l) gauged WZW model. Even so, this may also 
be misleading, as we will discuss in the final section. 



3. The type OA S-matrix 

The asymptotic states of type OA string theory are perturbations of the tachyon 
field, realized in the matrix description as disturbances of the eigenvalue density. To 
evaluate the S-matrix, one is instructed to perform a kind of LSZ reduction on eigenvalue 
density correlators, taking the point of evaluation of the corresponding fermion bilinears 
to |A| — > oo, t — > ±oo and extracting the coefficient of the asymptotic behavior. The 
connected correlation function of the asymptotic eigenvalue density perturbations consists 
of ring diagrams of the corresponding fermion bilinears; an example of 1 — > N scattering 
is depicted in figure 4. 
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Figure 4. 1 —>■ N scattering in the matrix model. The incoming tachyon (dashed 
line) splits into a particle (blue) and hole (red), which reflect off the potential wall. The 
particle-hole pair decays back toward the Fermi surface via the emission of final state 
tachyons. 



A recipe for the result of this calculation was given in |3^] (and reviewed in ||). One 
picks from the products of fermion propagators the terms in which the wildly oscillating 
phases exp[|A 2 ] cancel; one then extracts the coefficients of the asymptotic behavior of 
the density field 

p(X,t) ~ ~exp[-iw(t±log|A|)] (3.1) 

to obtain the scattering matrix element. For each fermion reflecting off the potential wall 
near the origin in A space, there is a factor of the reflection coefficient R(u) for particles 
and R*{u) for holes. The reflection coefficient of type OA is given by! 



R(u) 



h U 

4 16 P 



r(| + | + ^-^) 

r(| + | — iu + ifj,) 



(3.2) 



In comparing to worldsheet results, there are in addition certain "leg-pole factors" 
relating the matrix model density perturbation to the string tachyon []35|]. These are 
energy-dependent phases for on-shell scattering processes, which therefore affect quantities 
such as the time delay of an outgoing particle. However, they will not affect the outgoing 



We set a' = 2. 
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particle distribution as a function of energy, which is the object we study here. Therefore 
we henceforth ignore these factors .0 

Our operating hypothesis, as explained in the last section, is that any sufficiently 
energetic incoming state causes the formation of a black hole. We will thus concentrate on 
the 1 — > N amplitude for simplicity (a partial analysis of the 2 — > N amplitude is contained 
in the Appendix). The momenta of the external states then divide into two groups: The set 
S that attach to the propagating fermion line, and therefore have positive frequency; and 
the set S that attach to the propagating hole line,i and therefore have negative frequency. 
According to the prescription outlined above, the scattering amplitude is 

A N {u>\u>i)= Yl P (%R(u-Z)R*(-0 (3-3) 

Sc{l,-",iV} us 

where u s = J2i=i u ji and \ s \ = k for S = {j u ■ ■ -,j k }- 

Plugging the explicit form of reflection coefficient ( |3.2| ) into ( |3.3|) , An is rewritten as 



An 

where we defined F(£) by 



2 x — _ 

T"l6 +/1 



/ <%F(u-OF(0 (3.4) 



_ r(| + i + i|) 
m ~ r(| + I- i?) - (3 ' 5) 



,2 



Dropping the overall phase factor | ^ — + /j, \ IUJ and changing the integration variable 
as £ — > ^r^, the amplitude A n in ( |3.4j) becomes 



The leg-pole factors could potentially alter the result when considering a wave packet of 
incoming tachyon. The width Suj of the wavepacket determines the localization in position space. 
One can choose the width much smaller than the string mass, so that the leg pole factor is constant 
over the width of the wavepacket in frequency space, and yet much larger than the inverse of the 
time scale T that the entire wavepacket spends in the strong coupling region, behind the point in 
4> where the low energy field equations would predict that a horizon forms. The variation of the 
amplitude is smooth over the width of such a wavepacket, and so the result of scattering such a 
localized wavepacket will not differ substantially from that of a pure plane wave. 

5 There is always one tachyon emitted at the point where the particle and hole lines meet in 
the final state; in our conventions, this tachyon belongs to the set S. 
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In |3"2]1 , a saddle point approximation was applied directly to this expression (|3.6| ) of 
amplitude. However, it is more useful to perform a summation over S before doing the 
saddle point approximation. For this purpose, we utilize the trick of inserting "1" into the 
integral: 

1 = f°° du5(u - £) = r ^*e ft ("-€). (3.7) 



Then ( |3.5|) becomes 



pui-2n poo 



— g — 2/x J LOs~ UJ s~'^t 1 J — oo ^" 



dvdt i+(,,i_o,,\ „, s 1 



(3.J 



2tt w -it 



where we defined /(£) = Now we can perform the summation over S by 

using the relations ^2 S (— l)' 5 ' = and 

N 

^(_l)|S| e «(«s-«s) = JJ( e 4 *^ - e"^)- (3-9) 
s i=i 

Finally, we arrive at a compact form of A n 

/oo N 
^G(t)e 2lt ^ Jj2sin^t (3.10) 

where G(t) is given by 

Note that G(t) is an even-function of t. The parameter t can be interpreted physically as 
the time difference between the bounce of the particle and the bounce of the hole. 

In the following, we will set \x = for simplicity (its effects are easily restored). Our 
expression ( |3.10| ) reproduces the well-known selection rule that when \x = the amplitude 
vanishes for N = even [|36| |. Below we will study the high energy behavior of the non- 
vanishing N = odd amplitude: 



An = i 



f°° dt N 

N-i j jG{t) Y[2ajiu jt . (3.12) 
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3.1. Zeroth order approximation 

Suppose G(t) is independent of t: G(t) = G(0) Vt. Then the amplitude ( |3.12| ) becomes 

roo jx N 

A N = G^i' 1 / — TT(e^ - e~ lt ^) 
Jo t j=1 

f°° rU 

s Jo t (3.13) 



= |G(0)^(-l)l 5 le( 



= ttG(0) ^(-1)1^-2^) , 

s 

where e(x) is the sign of x and 9(x) is the step function. In the third equality of ( |3.13| ), 
we used the formula J °° * sma ^ = f e ( a )- This expression (|3.13| ) is the same as eq.(8.6) 
in [j32~l , which was obtained by applying the saddle point method directly to the integral 
( |3.6|) . As is obvious from the derivation of (|3.13| ), this approximation does not capture an 
essential part of reflection, contained in the t-dependence of G(t). We shall see shortly that 
G(t) is approximately constant, but only over a finite interval t < to ~ 0(uj~ 1 / 2 ). It turns 
out that the integral (|3.12| ) has most of its support in the region t > to when iV ^> y/u; a 
more detailed analysis is thus required. 



4. Saddle point approximation of reflection 

From the expression (|3.12j ), one can see that the amplitude vanishes linearly as u)\ — > 0, 
which is an example of the low energy theorem for soft tachyon scattering [|32| . Therefore, 
we expect that the dominant contribution is localized in the middle of the phase space away 
from the boundary {u>i = 0}, and when N and u> are sufficiently large the amplitude An 
is sharply peaked around the mean value (ui) = |4. This in turn implies that the saddle 
points of the t-integral are located at sm(uJi)t = ±1, i.e. t = t\ = ^L-)^ = 1, 2, • • We 
checked numerically that the peak around the higher saddle points t;>i easily disappears 
when UiS fluctuate around |4, and thus we conclude that the first saddle point t\ = \- 
gives the dominant contribution to the i-integral (|3.12| ). 

Numerical evaluation of G(t) shows that it is nearly constant for t < to ~ O(o; -1 / 2 ), 
and exponentially decaying for t > to (see figure 5 and below). This in turn leads to 
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a change in the characteristic behavior of the scattering amplitude for iV ^> y/uj and 
N <C y/u; we therefore analyze these two regimes separately. 



IGI 




10 11 12 13 14 



co 1/2 t 



Figure 5. Numerical evaluation of \G{t)\ for u = 400, q = 0. The curve for lu 
800, q — is virtually identical apart from a rescaling of the vertical scale by y/2. 



4-1- N ^> y/u case 

In the case N y/u, the saddle point t\ ~ — is much larger than the characteristic 
time scale to ~ u~^ of G(t). Therefore, the exponentially decaying behavior of G(t) for 
t 3> to is relavant for the evaluation of the t-integral ( |3.12| ). This behavior of G(t) can be 
reproduced by closing the contour of the v integral ( |3.11[ ) in the upper or lower half of the 
complex z/-plane (depending on the sign of £), and picking up the contribution from the 
poles of the V- functions in (|3.5|) . In this way G(t) is estimated as 



n=0 

2sinh7ru; 
iri 

2sinh7ru; 

ivi 



{-l) n r(l + q + iu + n) 
n\ T(l + q + n)T(-iu) — n) 



e -(i +g+fo , + 2n)|t| 1 r(l + g + w + n)T(l + iu + n) 

n=0 



-{X+q+iu)\t 



T(l + q + iu)T(l + iu) 



r(i + ?) 



2 Fi(l + q + iu, 1 + iu, 1 + q; e _2|t| ) . 

(4.1) 

We should emphasize that this computation of closing the contour does not give an exact 
answer of G(t) since the integrand of G(t) does not decay as v — > ±oo. 
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The n = term in the summation (|4.1|) gives 



r(i + g + w) (1+,+^ 
^ w r(i + ? )r(-zu;) e 



(4.2) 



This term comes from poles in ( |3.11|) a,t uj ± u = i(l + q) , where almost all the energy of 
the initial tachyon is carried by either the particle or by the hole. Including the higher 
order terms in the sum does not substantially alter the results. 

Now the t-integral in fl3.12| ) can be evaluated by expanding the integrand around the 
saddle point at t = ^ and Ui = jf? 



t 



TV 



N 



2u 



+ u, 



U>4 



N 



JV 

E 

i=i 



0, \u\, \eA <S 1 



(4.3) 



By performing the Gaussian integral over the fluctuation w, An is found to be 
An ~ 



r(l + q + iu) _ 2N 



T(l + q)T(-iu) 
T(l + q + iu) 

~ r(i + ?)r(-i W ; 

Here we assumed q <C u. 



du e -(l+9+^)(if +u) exp 



71- N 

2uj 



2N 



-u 



2 \2uo ) ^ 1 



2 N e -(l+q+iuj)^ 



N e~ N / 2 



-exp 



m 2 ^ 



(4.4) 



4.2. N <^ v /cJ Case 

In the opposite case 1 < JV < y/u, the relevant behavior of G(t) is the plateau region 
in t <C to (see figure 5). This behavior of G(t) can be obtained by expanding the zv-integral 
( |3. 11|) around v = 0. Physically, the point v = corresponds to the situation where both 
particle and hole carry the same amount of energy When v is small, the integrand in 
Q3.11|) is well approximated by a Taylor expansion in v 



Git) 



2n V2 



i{u) 2 + q" 



-F — 
2 



4uj L 



(4.5) 



with the saddle point 



u 2 + q 2 
2u 



(4.6) 
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This approximation is self-consistent if is smaller than the width of Gaussian in (4.5), 

(4.7) 



uj 2 +q* 



t < 



UJ 



UJ 



q 2 + uj 2 



Outside of this range, G(t) is highly suppressed due to the oscillatory factor e ltu . Therefore, 
G(t) is approximated as 



• ' 2 V Jrw 



. P (|)V^ 9(t0 _| t |) 



where 



to = 7 



Uj 



(4.f 



(4.9) 



q 2 + uj 2 

7 is a numerical coefficient independent of q and uj. Note that the time difference t is very 
small t < uj~^ in the high energy limit. Putting it all together, in this regime the high 
energy 1 — ► N amplitude is given by 



A N = -r 



tvuj V2 



2 f to dt 
— ( 
t 



N 



Jj2sinwji. (4.10) 



This approximation is consistent when ^ -C beyond this, the Taylor expansion ( |4.5| ) 
breaks down and it is better to use ( |4.1|) . Assuming q <C u, An becomes (up to a phase 
factor) 



dt 



N 



An ~ v - / — e ' 1 



j^J 2 sinwjt 



i=i 




A 



exp 



2 V 2^ / ^ 



(4.11) 



Here we neglected the it-dependence coming from the factor e , since it is subleading 
when - < 1. 



5. The number distribution at high energy 

One of the physically interesting quantities we can calculate is the number distribution 
Vn of final state tachyons in the 1 — ► N scattering, which is defined by 



1 1 



TV 



duji 



1=1 I 



UJ; — UJ 



AN{w\oJi 



(5.1) 
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Vn for the low energy and weak coupling scattering has been studied in [B2| and it was 
found that Vn is almost a Poisson distribution. We are interested in the high energy 
behavior of Vn- 

If we substitute the naive step function approximation of A n (|3.13|) into the definition 
of Vn (pTT|) , it turns out that the resulting c^-integral can be performed exactly: 

(2n) 2k - 1 (2 2k -l)\B 2 k\ , \ 

p — ~ mz (5 ' 2) 

Here we used the approximation |7rG(0)| ~ a/iF as m Note that ( |57^ ) is a mono- 

tonically increasing function of N = 2k — 1. It is argued []32| that the approximation 
( p.!3| ) breaks down near the boundary of step functions. However, one can argue that the 
contribution from the dangerous boundary region is negligible for small N, so we expect 
( |5.2|) to be a good approximation for sufficiently small N. 

When N becomes large, ( |5.2[ ) cannot be trusted anymore and we should use the 
approximation (fOj) or ( |4.11|) , depending on the regime iV 3> y/u or 1 C < y/u. In 
both cases, the Wi-mtegral can be approximated by the Gaussian integral around the saddle 
point (ui) = % 



exp 



ttN \ \ - --, 
^7 I> 



N 

Therefore, in the regime iV 3> y/u, Vn is given by 

1_ 

N\ ' 

and similarly in the regime iV -C y/uj 



2 \ N ~ 1 N 



71/ UJ 2 



(5.3) 



N 

ttN 



V «~W.W{^ ] e_ ™" (5 ' 4) 



^-^-.(Vi" ■ (5-5) 



The discrepancy in these two formulae for N ~ y/u is sub leading in N, and can be 
attributed to the different approximations made in the two regimes. 



Curiously, there is a nice answer for the generating function of this result: 

VkX k = (2a;) -1 tan(7rx) . 



fe=i 



We also find it remarkable to discover yet another appearance of the seemingly ubiquitous Bernoulli 
numbers. 
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6. Discussion 

The main observation is that the probability to have iV outgoing particles is a mono- 
tonically and steeply decreasing function of N, at least when N is sufficiently large. On 
the other hand, ( |5.2| ) suggests the probability grows with N for iV small; indeed, in the 



probability to have more than one particle in the final state at high energy was estimated 
by a unitarity argumentB to be (noting that V\ ~ according to ( |5.2| )) 



V N>1 = 1 - Vi ~ 1 - . (6.1) 

Taken together, our results suggest that the distribution is peaked for some reasonably 
small number of outgoing particles, larger than one. Note that equation ( |5.2| ) gives Vn ~ 
0(1) for N ~ logo;; we therefore expect Vn to be peaked in this regime. 

There is no characteristic feature for N ~ uj that might indicate the presence of a black 
hole or similar object as an intermediate state. Instead, the configuration that dominates 
high energy scattering is the motion of a single particle-hole pair bouncing off the potential. 
For N ^> y/uj, the energy is carried predominantly by either the particle or the hole; for 
N <C y/u, it is shared roughly equally by each. 

A qualitative explanation for the rapid decrease of Vn with iV comes from the fact 
that the iV tachyon state is built from fermion bilinears The initial state consists of 

one such bilinear, which then picks up a phase from reflection off the potential ( |1.4j ) . The 
final state is a product of N such bilinears, and the overlap with the evolved initial state 
picks out the terms in ip *ip • • • ip where all but two of the fermion operators act on one 
another. This is a relatively small fraction of the overall possibilities, since this operator 
can create up to iV particle-hole pairs when acting on the vacuum. The initial state selects 
only the part with one particle- hole pair. 

6.1. Comments on conservation laws 

The free fermionic character of the matrix model leads to an infinite set of conserved 
quantities; the time-independent ones are 

Qt= f dee l b\b e (6.2) 



7 Unitarity also tells us that the probability to send in a single initial tachyon and then find N 
tachyons in the final state, is the same as that for the time reversal prosess of iV — > 1 scattering. 
Intuitively, it seems unlikely to have a single tachyon state as the final state of N tachyon scattering 
with large N, unless we fine-tune the initial state. 
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where 6, are fermion creation and annihilation operators. Typically, black holes have 
no hair, suggesting that their conserved charges are Qi = u, Qe = O.i States with these 
properties in the free fermion Hilbert space consist of a macroscopic number of very soft 
tachyon excitations - more or less a coherent state of soft tachyonsJ! On the other hand, 
our one-particle initial state has Qi ~ u/. In order to form a black hole as an intermediate 
state, an initial state consisting of a small number of high energy tachyons would have to 
shed these higher conserved charges, e.g. by radiating a shell of outgoing tachyons carrying 
them away, leaving behind an intermediate state with the characteristics of the black hole. 

The phase space available to soft particle-hole pairs is not large, and the only objects 
with the charges Qi ~ oj5i^ are collections of soft fermion excitations of energies {wj}, 
i = 1, N, relative to the filled Fermi sea;0 these yield the conserved charges 



N N 

Q 1 =^±u i =u ■ Qi = Y,±ui~o{N 1 - i ) , £>1 (6.3) 

i=l i=l 

with the sign + for particles and — for holes. Since there is no exponentially large phase 
space of black hole states, the evolution is not drawn towards the formation of an in- 
termediate state with the conserved charges associated to the black hole - there are not 
sufficiently many states to make visiting such a configuration likely. 

Now, to create an intermediate state corresponding to a finite deformation of the 
Fermi sea, such as a blob of eigenvalues separated away from the sea, one must consider a 
more general N' — * N scattering with N', N 3> 1, instead of the 1 — > N scattering studied 
in this paper. This is because the number of particle-hole pairs appearing in the N' — > N 
scattering is bounded from above by minjiV', N}. In particular, only a single particle-hole 
pair contributes in the 1 —>■ N scattering as we saw above, and therefore the initial state 
we considered has rather small overlap with the states carrying the charges expected for a 
black hole state. But if it happens that most high energy states do not scatter into black 
holes, then this casts doubt on the existence of any object that one should call a black 
hole. 

The intuition gained from our study of 1 — > N scattering shows why this soft-excitation 
state will not have the properties of a black hole. The Hawking radiation of a black hole 



Results to this effect were presented by A. Sen at Strings2004 [37]. See also [38]. 
9 This feature is in accordance with higher dimensional examples, where to form a macroscopic 
black hole requires the excitation of a macroscopic number of branes. 
10 So that uji < for holes and uJi > for particles. 
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consists predominantly of quanta with energy of order the limiting temperature T H = 
77- \ -=r. Such an outgoing state would have to be constructed from an intermediate state 
consisting of iV ^> u quanta, in order that the charges ( |673|) approximate those of a black 
hole, Qi = udifi'i again, the probability that N/u = n ^> 1 soft particle-hole excitations 
will make a Hawking quantum behaves like the overlap of a state created by (ip*i(j) n and 
one created by ip'if), which we expect to be a rapidly decreasing function of n. 

6.2. Whither the 2D black hole? 

The evidence seems to suggest that there is no such object as the 2D black hole in 
string theory. The question then is why we were misled into thinking so by the seemingly 
exact coset construction of a worldsheet CFT describing strings propagating in a black 
hole background.^ Consider the Euclidean continuation of (|1.3|) , in which the putative 
Euclidean black hole geometry is a capped semi-infinite cylinder often referred to as the 
"cigar" . This geometry is strongly curved in the vicinity of the tip of the cigar for the regime 
of interest, k ~ 2. In fact, there is a conjectured strong-weak coupling duality |H| , fi2"| , |39| 



between the 2D black hole and Sine-Liouville theory, i.e. the tachyon background 

T = \e^cos9 (6.4) 



describing a condensate of vortices, where 9 is the (axial) target Euclidean time coordinate, 
and Q = (k - 2)" 1 / 2 . The Euclidean cigar geometry is a weakly coupled sigma model for 
Q — > 0, while the Sine-Liouville background is weakly coupled for Q — > 00. The regime of 
interest for 2D black holes is the strongly coupled one, where the Sine-Liouville description 
is more appropriate, and one cannot use the cigar geometry to extract reliable predictions. 
This observation is consistent with the fact that the Hagedorn thermodyamics of high mass 
black holes is not to be found in the matrix model. The Sine-Liouville background appears 
to describe instead a classical tachyon background of the 2D string, consistent with the 
picture deduced from the consideration of conserved charges in the previous subsection. 

Perhaps the proper lesson to be drawn from all the known facts is that high energy 
scattering is dominated by brane physics. In higher dimensions, brane physics is black 
hole physics - the degrees of freedom on brane intersections lead to a large degeneracy 
of black hole intermediate states, which dominate the S-matrix. In two dimensions, the 
high energy scattering is still dominated by branes - the DO-branes that are the eigenvalue 
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We thank David Kutasov for the following argument; see [39,4C] . 
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fermions of the matrix model - but the kinematics is so restrictive that they do not have 
a black hole density of states. Instead of forming a black hole, a highly energetic tachyon 
in a sense becomes a DO-brane that bounces high off the background potential of figure 2, 
and then fragments into a small number of very energetic tachyons in the final state. 

Acknowledgments: We wish to thank David Kutasov and Jongwon Park for useful 
discussions, and Sumit Das and Sergei Gukov for correspondence. This work is supported 
in part by DOE grant DE-FG02-90ER40560. 



Appendix A. The 2 — > N amplitude 



The 2 — > N amplitude is also written in terms of the function G(t) in ( |3.11| ). When 
H = and N = even, the 2 — > N amplitude is given by 



A 



2^N 



UJ UJ 

2' 2 



1 



7T 



dt 



ds 



G u (t) 



N 



N 



j^J 2i sin ui (s + t) — TT 2i sin uji (s — t) 



i=i 



AUB={1, — N} 



Y[ 2i sin oji (s + t) — J^J 2i sin ui (s — t) 



lie a 



leA 



n 2i ~ x si 

jeB 



sin ujj s 



Here f * g(t) is the convolution 

f*g(t) 



duf(u)g(t - u) . 



(A.1) 



(A.2) 



Qualitatively, this amplitude will have a similar structure to the 1 — > N amplitude. Both 
G and G * G only have appreciable support for t < (^(a; -1 / 2 ), and the sine factors are 
qualitatively similar to the 1 — > N result. Thus we expect qualitatively similar behavior, 
namely a monotonically decreasing probability for finding N particles in the final state as 
a function of N. 
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